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In  this  study,  the  dynamic  detonation  parameters,  namely,  the  critical  tube  diameter  and  critical  energy 
for  direct  initiation,  of  C2H2/N20/Ar  mixtures  were  measured  at  various  initial  conditions.  Using  chemical 
kinetics  with  the  Konnov  mechanism  and  experimental  measured  data,  a  simple  correlation  to  evaluate 
the  critical  tube  diameter  of  C2H2/N20/Ar  detonation  is  developed,  given  by  dc  =  594. 8<p0623 
(1  -XAr)0  2176(p/po)“0  0246zl/  which  is  applicable  for  stoichiometric  mixtures  with  initial  pressures  ranging 
from  50  to  130  kPa,  equivalence  ratios  from  0.625  to  2.5  and  stoichiometric  mixtures  with  maximum  per¬ 
centage  of  argon  dilution  up  to  50%  at  the  atmospheric  pressure.  By  combining  this  correlation  function 
with  a  theoretical  model  based  on  a  simple  work  done  concept,  the  critical  energy  for  direct  blast  initi¬ 
ation  can  be  predicted;  and  the  theoretical  predictions  are  found  to  be  in  good  agreement  with  the  critical 
energy  measured  experimentally.  To  assess  the  detonation  sensitivity,  the  critical  energy  results  for 
direct  initiation  of  detonation  in  C2H2/N20/Ar  mixtures  are  compared  with  those  in  H2/N20/Ar  mixtures 
[Zhang  et  al.,  Int.  J.  Hydrogen  Energy  39  (2011)  5707-5716].  The  results  indicate  that  C2H2/N20/Ar  is 
more  sensitive  than  H2/N20/Ar  at  the  same  initial  conditions.  This  is  also  supported  qualitatively  by 
the  chemical  kinetic  calculation  which  shows  that  the  ZND  induction  length  scale  for  C2H2/N20/Ar  is  rel¬ 
atively  smaller  than  that  of  H2/N20/Ar  mixture.  Lastly,  the  effect  of  the  oxidizer  on  the  detonation  sensi¬ 
tivity  is  studied  by  comparing  the  critical  energy  between  C2H2/N20/Ar  and  C2H2/02/Ar  mixture.  It  is 
found  that  a  higher  energy  is  required  to  successfully  initiate  a  spherical  detonation  in  C2H2/N20/Ar  than 
in  C2H2/02/Ar  mixtures  and  equivalently,  the  results  agree  qualitatively  with  the  ZND  induction  length 
analysis. 

©  2012  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Nitrous  oxide  (N20)  is  a  stable  compound  that  is  comparatively 
unreactive  to  many  substances  at  room  temperature.  Its  use  has 
been  considered  in  medical  application  for  its  anesthetic  effects 
and  rocket  propulsion  research  as  a  cold  gas  propellant  or  oxidant 
[1,2].  When  heated  sufficiently,  however,  N20  is  an  energetic 
oxidizer;  it  first  decomposes  endothermically  to  N2  and  O  by  the 
initial  decomposition  step:  N20  (+M)  <-►  N2  +  O  (+M)  and  subse¬ 
quently  energy  release  can  be  driven  by  oxygen  atoms  combina¬ 
tion:  0  +  0  (+M)  <-►  02  (+M).  Thus,  N20  mixed  with  other  gaseous 
fuels  can  potentially  lead  to  various  accidental  explosions  [3,4]. 
Therefore,  it  is  necessary  to  understand  the  subject  of  nitrous  oxide 
explosion  properties  and  to  assess  its  potential  hazards. 

For  the  past  decades,  many  studies  have  been  devoted  to  ex¬ 
plore  the  combustion  properties  such  as  ignition  delay  times,  flame 
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burning  speeds,  flammability  limits  and  detonation  cell  sizes  in 
hydrogen-nitrous  oxide-diluent  mixtures  [5-13].  These  studies 
were  mostly  motivated  by  the  concern  of  hazards  posed  by  the 
accidental  explosions  in  these  H2/N20  gaseous  mixtures  that  are 
continuously  generated  from  stored  wastes  of  military  nuclear  pro¬ 
cesses  [14-16].  Unlike  common  hydrocarbon  fuel-oxygen  or 
hydrocarbon  fuel-air  mixtures,  the  detonation  phenomenon  in 
hydrocarbon  fuel-nitrous  oxide  mixtures,  e.g.,  C2H2/N20  consid¬ 
ered  in  this  study,  has  not  been  widely  investigated.  Such  studies 
could  have  an  impact  of  improving  the  safety  in  industrial  sectors 
that  deals  with  nitrous  oxide  from  waste  treatments,  emissions 
from  energy  sources  or  chemical  processes  [17],  and  the  applica¬ 
tion  of  nitrous  oxide  as  an  oxidizer,  for  example  in  chemical  rocket 
engines  [2].  Among  different  detonation  parameters,  evaluation  of 
critical  diameter  and  critical  energy  for  direct  initiation  provides 
direct  means  to  assess  the  detonation  sensitivity  of  a  given  mixture 
[18,19]  and  yet,  data  on  these  dynamics  parameters  for  hydrocar¬ 
bon  fuel-nitrous  oxide  mixtures  are  rather  scarce.  Although  deto¬ 
nation  cell  size  is  also  often  considered  as  an  important  dynamic 
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parameter  to  quantify  the  detonation  sensitive  of  a  given  explosive, 
it  was  found  that  fuel-N20  mixtures  tend  to  be  very  unstable  with 
highly  irregular  cellular  structures.  For  example  in  C2H2/N20/Ar 
mixtures,  even  with  the  argon  dilution  up  to  50%,  Libouton  et  al. 
[20,21]  observed  “substructure”  with  cell  sizes  of  multiple  scales, 
which  makes  the  determination  of  the  cell  size  for  C2H2/N20/Ar 
mixtures  very  difficult. 

In  this  study,  we  focus  on  the  C2H2/N20/Ar  systems  as  an  exam¬ 
ple  in  investigating  the  effect  of  nitrous  oxide  on  the  sensitivity  of 
detonation  formation  in  hydrocarbon  fuel  mixtures.  Critical  tube 
diameter  and  critical  energy  for  direct  initiation  data  for  C2H 2\ 
N20/Ar  mixtures  are  experimentally  obtained  at  different  initial 
conditions  to  look  at  the  effects  of  the  initial  pressure,  equivalence 
ratio  and  the  amount  of  argon  dilution  on  the  detonation  dynamic 
parameters.  These  results  are  then  compared  with  those  for  H 2\ 
N20/Ar  to  study  the  influence  of  hydrogen  as  a  fuel,  and  compared 
with  C2H2/02/Ar  mixtures  to  observe  the  detonation  sensitivity  if 


nitrous  oxide  is  being  used  as  an  oxidizer.  To  investigate  the  deto¬ 
nation  sensitivity  from  chemical  kinetics  consideration,  an  induc¬ 
tion  zone  length  analysis  is  also  carried  out  in  parallel  to  support 
the  experimental  observations. 

2.  Experimental  details  and  chemical  kinetic  analysis 

2.1.  Experiment  setup 

Measurement  of  the  critical  tube  diameter  and  critical  energy 
were  both  carried  out  in  a  modified  spherical  bomb  apparatus  as 
shown  in  Figs.  1  and  3,  respectively.  The  original  spherical  chamber 
is  20.3  cm  in  diameter  and  5.1  cm  in  wall  thickness.  In  both  exper¬ 
iments,  each  explosive  mixture  was  first  prepared  beforehand  in  a 
separate  vessel  by  the  method  of  partial  pressures.  The  storage 
pressure  was  kept  at  200  kPa  and  the  gases  were  allowed  to  mix 
in  the  vessel  for  at  least  24  h  in  order  to  ensure  homogeneity  for 
each  experiment. 

For  the  critical  tube  diameter  measurement  of  C2H2/N20/Ar 
mixtures,  experiments  were  carried  out  in  a  setup  consisted  of  a 
vertical  circular  steel  tube  connected  to  a  spherical  bomb  chamber 
(see  Fig.  1).  The  tube  diameters  were  varied  via  inserting  smaller 
diameter  tubes.  The  tube  diameters  used  in  the  experiment  were 
d  =  19.1  mm,  15.5  mm,  12.7  mm  and  9.13  mm.  The  initial  condi¬ 
tions  of  the  critical  tube  diameter  experiment  are  shown  in  Table  1. 
For  each  test,  the  detonation  was  initiated  by  a  high  voltage  spark 
ignition  source  at  the  top  of  the  vertical  circular  steel  tube  and  the 
wave  subsequently  transmitted  into  the  spherical  chamber.  A 
photo  probe  and  a  shock  pin  were  mounted  at  the  top  and  bottom 
of  the  spherical  bomb,  which  were  used  to  determine  the  time-of- 
arrival  (TOA)  signal  of  the  wave.  The  successful  or  unsuccessful 
detonation  when  it  transmits  into  spherical  chamber  can  be  esti¬ 
mated  by  the  TOA  signal  (see  [22]  for  further  details).  For  com¬ 
pleteness,  typical  traces  for  a  successful  detonation  transmission 
from  the  small  circular  steel  tube  (d  =  19.1  mm)  to  the  spherical 
chamber  with  relatively  much  larger  space,  i.e.,  for  d  >  dc,  in  a  stoi¬ 
chiometric  C2H2/N20  mixture  at  an  initial  pressure  of  p0  =  58  kPa  is 
shown  in  Fig.  2.  It  can  be  seen  from  Fig.  2  that  at  this  initial  pres¬ 
sure,  the  arrival  time  of  the  wave  at  the  photo  probe  is  206.8  ps 
and  332.1  ps  at  the  shock  pin.  The  velocity  of  the  wave  is  calcu¬ 
lated  to  be  2021.3  m/s  and  2002.4  m/s  in  the  vertical  tube  and 


Voltage  vs  Time 


Fig.  2.  Arrival  time  trace  showing  a  planar  detonation  successfully  emerged  from  the  small  steel  tube  d  =  19.1  mm  into  the  large  spherical  chamber  for  the  stoichiometric 
C2H2-5N20  mixture  at  an  initial  pressure  of  58  kPa. 
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Table  2 

Initial  conditions  used  in  the  direct  initiation  experiment. 


Mixtures 

Initial  pressure 
(kPa) 

Equivalence 
ratio  (q>) 

Argon  dilution 

(%) 

C2H2/N20/Ar 

50-200 

0.56-3.0 

0-50 

H2/N20/Ar 

70-200 

0.5-1. 2 

0-30 

C2H2/02/Ar 

5-25 

0.625-4 

0-70 

Slender  Coaxial  Electrode 


PCB(113B24) 


Heating  Tape 

Fig.  3.  Experimental  setup  for  the  direct  blast  initiation  experiment. 


Table  1 

Initial  conditions  used  in  the  critical  tube  diameter  experiment  for  C2H2/N20  mixture. 


Equivalence  ratio  (cp) 

Initial  pressure  (kPa) 

Argon  dilution  (%) 

0.625-2.5 

38-200 

0-50 

spherical  chamber,  which  is  approximately  90.8%  and  90%  of  the  CJ 
detonation  velocity,  respectively.  It  shows  that  at  an  initial  pres¬ 
sure  of  58  kPa,  the  tube  diameter  is  above  the  critical  value,  thus 
the  planar  detonation  can  successfully  transmit  into  a  spherical 
detonation.  On  the  other  hand,  for  an  unsuccessful  transmission 
with  an  initial  pressure  lower  than  the  critical  value,  e.g., 
p0  =  52  kPa,  although  a  detonation  wave  still  formed  in  the  vertical 
tube,  the  detonation  failed  after  exiting  into  the  free  space.  At  this 
initial  pressure,  the  wave  first  reached  the  photo  probe  at 
t  =  211.2  jis,  and  the  shock  pin  at  t=1180ps.  With  these  TOAs 
from  the  traces,  the  velocity  of  wave  in  the  vertical  tube  is  roughly 
about  90%  CJ  and  in  the  spherical  bomb  the  velocity  is  estimated  to 
be  around  11.4%  CJ,  which  clearly  indicates  failure  of  the  wave.  For 
each  successful  and  unsuccessful  initiation  of  spherical  detonation 
at  least  3  shots  were  repeated  to  confirm  the  critical  pressure  that 
can  form  a  spherical  detonation  at  each  tube  diameter. 

As  for  the  critical  energy  measurement,  direct  initiation  of 
spherical  detonations  was  achieved  via  a  spark  discharge  from  a 
high  voltage  and  low  inductance  capacitor  inside  the  closed  spher¬ 
ical  bomb  chamber  shown  in  Fig.  3.  In  brief,  the  ignition  system 
consists  of  a  high  voltage  power  supply,  capacitor  bank,  a  gap- 
switch,  a  trigger  module  (TM-11A)  and  a  slender  coaxial  electrode 
mounted  on  top  of  the  explosion  sphere.  At  the  end  of  this  slender 
electrode  there  is  a  3.5  mm  spark  gap  through  which  the  energy  is 
delivered  inside  the  chamber  through  the  ignition  circuit.  A  PCB 
piezoelectric  pressure  transducer  is  mounted  in  the  wall  of  the 
chamber  to  measure  the  time  of  arrival  of  the  wave  front.  From 
the  arrival  time  of  pressure  signal  it  can  be  determined  whether 
there  is  a  successful  direct  initiation  of  detonation  or  not.  The  pro¬ 
cedure  to  distinguish  detonation  initiation  as  well  as  details  to  esti¬ 
mate  the  actual  spark  discharge  energy  from  the  ignition  system 
can  be  found  in  authors’  previous  studies  [23-28].  In  this  study, 
a  wide  range  of  initial  conditions  of  C2H2/N20/Ar,  H2/N20/Ar  and 
C2H2/02/Ar  mixtures  were  used,  which  are  summarized  in  Table  2. 

2.2.  Chemical  kinetics  modeling 

To  investigate  the  detonation  sensitivity  of  mixtures  under  dif¬ 
ferent  conditions,  some  aspects  of  chemical  kinetics  (i.e.,  induction 


length  analysis)  was  also  considered.  The  induction  length  analysis 
often  provides  some  insights  to  explain  qualitatively  the  variation 
of  critical  initiation  energy  and  critical  tube  diameter  observed  in 
the  experiments.  Such  analysis  was  commonly  used  by  a  number 
of  detonation  studies  in  literature  [29-39]  where  the  one-dimen¬ 
sion  ideal  ZND  model  is  used  to  calculate  the  characteristic  induc¬ 
tion  length  scale.  The  ZND  induction  length  is  generally  referred  to 
the  length  of  the  thermally  neutral  period,  computationally  de¬ 
fined  by  the  distance  from  the  leading  shock  to  the  location  of 
the  maximum  temperature  gradient  in  the  ZND  detonation 
structure. 

It  is  however  important  to  note  that  the  ZND  analysis  using 
chemical  kinetics  does  not  take  into  account  the  dynamic  effect, 
for  which  detonation  cell  size  2  should  represent  a  more  suitable 
characteristic  length  scale  to  describe  the  dynamic  detonation 
structure.  In  fact,  all  the  previous  detailed  ZND  analyses  show  that 
there  appears  no  direct  universal  proportionality  between  dynamic 
parameters  with  the  ZND  reaction  length  scale  or  equivalently  cell 
size  with  induction  zone  length  [30-37].  Often  the  variation  of  the 
proportionality  coefficient  must  be  suitably  determined  using 
some  correlation  models  [38,39]  in  order  to  carry  any  further 
quantitative  analysis  on  different  effects  of  mixture  composition 
and  initial  conditions.  In  any  case,  however,  dimensional  consider¬ 
ation  must  suggest  that  the  critical  tube  dc,  cell  size  2  and  ZND 
induction  length  A  should  be  scaled  with  each  other,  e.g.,  2  ~  A  •  A, 
although  the  proportionality  coefficient  A  may  not  be  constant 
depending  on  different  mixture  compositions  and  initial  condi¬ 
tions.  For  the  C2H2/N20/Ar  mixtures  investigated  in  this  work, 
the  variation  of  A  is  nevertheless  found  to  be  within  the  same  order 
of  magnitude.  (This  is  supported  by  the  correlation  presented  in 
the  next  section  (see  Eq.  (1))  from  which  one  can  show  that  the 
variation  is  only  substantial  near  the  rich  composition  limit, 
(p  >  2.0.)  The  same  observation  applies  for  the  H2/N20/Ar  using  a 
similar  correlation  given  in  [25]  for  detonation  cell  sizes.  The  var¬ 
iation  of  the  corresponding  proportionality  coefficient  A  therefore 
should  not  affect  the  qualitative  trends  of  the  results  described 
by  the  ZND  analysis  in  this  work  and  the  chemical  kinetic  calcula¬ 
tions  should  be  useful  to  elucidate  and  as  well  capture  qualitatively 
the  effects  of  mixture  composition  and  initial  pressure  on  the  crit¬ 
ical  tube  diameter  and  the  critical  initial  energy. 

All  the  chemical  kinetic  computations  performed  in  this  study 
used  the  CHEMKIN  package  [40].  For  C2H2/N20/Ar  and  C2H2/02/ 
Ar  mixtures,  the  Konnov  chemical  kinetic  mechanism  [41]  was 
used.  The  version  of  Konnov  mechanism  used  in  this  work  contains 
1016  reactions  and  122  species  (Release  0.3).  This  version  of  the 
Konnov  mechanism  has  been  validated  with  experimental  data 
for  oxidation,  ignition,  and  flame  structure  of  small  hydrocarbon 
and  also  cases  with  nitrous  oxide  [41].  It  was  largely  revised  using 
the  data  from  [42]  for  acetylene  oxidation.  In  addition,  other  ver¬ 
sion  of  Konnov  mechanism  is  often  chosen  for  detonation  analysis 
in  the  literature  [27,43-45].  Release  0.4  of  the  Konnov  mechanism 
was  also  validated  through  the  comparison  of  constant  volume 
explosion  simulation  and  available  shock  tube  ignition  data  for 
ethylene,  propane  and  hydrogen  [46].  For  mixtures  with  N20,  the 
sub-mechanism  of  nitrogen  chemistry  in  the  Konnov  mechanism 
is  developed  on  the  basis  of  the  widely  used  Miller-Bowman 
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mechanism  [47]  with  additional  reactions  and  updates.  It  is  worth 
noting  that  the  Miller-Bowman  mechanism  is  also  considered  in 
[6,13]  as  the  basis  and  modified  for  studying  detonations  in  mix¬ 
tures  containing  nitrous  oxide.  For  H2/N20/Ar  mixtures,  the  reac¬ 
tion  mechanism  constructed  by  Mevel  et  al.  [8]  was  used,  which 
had  been  well  validated  with  several  fundamental  combustion 
problems  including  ignition  delay  time  and  laminar  flame  speed 
[8,9],  and  prediction  of  detonation  properties  such  as  detonation 
cell  size  and  critical  initiation  energies  [10,25].  It  is  worth  noting 
that  this  chemical  kinetic  model,  constituted  of  203  reactions 
and  32  species,  has  also  been  constructed  from  the  mechanisms 
of  Konnov  [41]  and  of  Mueller  et  al.  [48].  Comparison  for  H2/N20 
mixtures  was  made  in  this  study  between  the  Mevel  et  al.  mecha¬ 
nism  and  the  Konnov  mechanism;  results  can  be  found  in  the  Sup¬ 
plemental  materials  of  this  paper.  It  is  found  that  the  difference  is 
not  substantial  and  qualitatively  the  variation  remains  essentially 
the  same.  Results  for  C2H2/02  mixtures  were  also  compared  be¬ 
tween  several  other  mechanisms  with  nitrogen  chemistry  with 
the  Konnov  model  and  the  difference  seems  reasonable,  although 
the  Konnov  mechanism  appears  to  estimate  a  smaller  ZND  induc¬ 
tion  zone  length.  Furthermore,  results  obtained  for  C2H2/02  mix¬ 
tures  using  different  versions  of  Konnov  mechanism  (Release  0.3, 
0.4  and  0.5)  are  very  close  to  each  other.  However,  for  C2H2/N20 
mixtures  there  seems  to  have  no  convergence  in  the  results.  Unfor¬ 
tunately,  to  the  author’s  knowledge,  there  is  no  experimental  data 
such  as  shock  tube  ignition  measurement  available  for  C2H2/N20/ 
Ar  mixtures  for  direct  comparison.  It  is  thus  difficult  to  identify 
which  chemical  kinetic  mechanism  is  particularly  suitable  for 
C2H2/N20/Ar  mixtures  and  used  for  this  study.  Therefore,  due  to 
its  reasonable  prediction  for  H2/N20/Ar  and  C2H2/02/Ar  mixtures 
and  for  consistency  purpose,  the  Konnov  mechanism  (Release 
0.3)  was  kept  for  all  the  calculations  in  this  work  for  the  acety¬ 
lene-based  mixtures  and  Mevel  et  al.  mechanism  for  H2/N20/Ar 
mixtures. 

3.  Results  and  discussion 


successfully  transmit  into  an  unconfined  spherical  detonation  - 
for  non-diluted  C2H2/N20  mixtures  at  different  composition.  The 
critical  condition  was  determined  by  carrying  out  the  experiment 
several  times  at  the  same  condition  (e.g.,  same  p0).  It  occurred  that 
at  initial  pressures  near  the  critical  value  for  a  given  tube  diameter, 
mixed  results  were  obtained  (i.e.,  sometimes  GO,  sometimes  NOGO 
at  the  same  initial  pressure  p0).  Therefore,  the  error  bars  shown  in 
the  plots  covered  the  range  of  the  mixed  results.  The  upper  and 
lower  error  bars  show  the  initial  pressures  above  or  below  which 
no  failure  or  successful  transmission  were  observed,  respectively 
within  the  number  of  shots  performed  at  each  condition.  From 
Fig.  4,  it  is  of  interest  to  note  that  C2H2/N20  at  (p  =  2.5  is  the  most 
sensitive  mixture  among  all  the  compositions  considered  here. 
Similarly,  Fig.  5  shows  the  experimentally  measured  critical  tube 
diameter  as  a  function  of  the  amount  of  argon  dilution.  For  50%  ar¬ 
gon  diluted  mixtures,  part  of  the  experimental  results  were  ex¬ 
tracted  from  Laberge  et  al.  [49],  which  follows  well  the  trend  of 
the  present  data  with  tube  diameter  smaller  than  20  mm.  In  gen¬ 
eral,  the  results  show  that  with  increasing  amount  of  argon  dilu¬ 
tion,  the  critical  tube  diameter  increases  consequently  at  the 
same  initial  pressure. 

In  these  plots,  correlations  are  also  included  for  qualitative 
comparison  and  to  further  elucidate  the  trend  of  variation.  The  cor¬ 
relation  is  based  on  the  consideration  that  the  detonation  front  can 
be  characterized  by  a  characteristic  reaction  zone  length,  and  the 
diffraction  process  is  governed  physically  by  the  tube  diameter. 
Dimensional  consideration  thus  suggests  a  correlation  between 
these  two  characteristic  length  scales.  Using  the  Konnov  chemical 
reaction  mechanism  for  C2H2/N20/Ar  mixtures,  the  induction  zone 
length  of  a  steady  ZND  detonation  structure  can  be  computed, 
which  can  be  used  to  scale  the  critical  tube  diameter  with 
dc  =  A  zl/.  Using  the  experimentally  measured  critical  tube  data, 
an  empirical  expression  for  computing  the  proportionality  param¬ 
eter  A  has  been  determined  -  a  mathematical  form  similar  to  those 
used  in  [10,25]  for  detonation  cell  size  correlation.  A  is  expressed 
as  a  function  of  equivalence  ratio  (p,  argon  diluent  mole  fraction 
XAr,  and  ratio  of  initial  pressure  to  standard  pressure  p/p°: 


3.1.  C2H2/N20/Ar  critical  tube  diameter  measurement  and  prediction 

Experimental  data  on  the  critical  tube  diameter  for  C2H2/N20 
mixtures  at  different  initial  conditions  are  presented  in  Figs.  4 
and  5.  Figure  4  first  shows  the  variation  of  tube  diameter  versus 
the  critical  pressure  -  above  which  the  planar  detonation  can 
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A  =  594.8(p°'623(l  -XAr)02176(p/pT0-0246  (1) 

For  the  above  correlation,  the  coefficient  of  determination  R2, 
and  the  maximum  deviation  between  the  correlated  values  and 
the  experimental  data  are  0.909  and  10.67%,  respectively.  Thus,  it 
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Fig.  4.  The  variation  of  tube  diameter  versus  initial  pressure  at  various  composi¬ 
tions  for  C2H2/N20  mixtures. 


Fig.  5.  The  variation  of  tube  diameter  versus  initial  pressure  at  different  amount  of 
argon  dilution  for  C2H2/N20  mixtures. 
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Fig.  6.  Critical  energy  as  a  function  of  initial  pressure  from  experimental  measure-  Fig.  7.  Critical  energy  as  a  function  of  equivalence  ratio  from  experimental 

ment  and  theoretical  prediction  for  stoichiometric  C2H2/N20  mixture.  measurement  and  theoretical  prediction  for  C2H2/N20  mixture  at  p0  =  100  kPa. 


can  be  seen  that  the  correlation  from  the  dc  =  A-AI  relationship  is 
in  fairly  good  agreement  with  the  experimental  data.  This  expres¬ 
sion  provides  a  reasonable  estimation  of  critical  tube  diameter  for 
C2H2/N20/Ar  mixtures  within  the  range  of  initial  pressures  from  50 
to  130  kPa  at  stoichiometric  condition,  equivalence  ratio  from 
0.625  to  2.5  and  with  maximum  percentage  of  argon  dilution  in 
stoichiometric  C2H2/N20/Ar  mixtures  up  to  50%  at  atmospheric 
pressure.  By  using  the  correlation  the  critical  tube  diameter  of 
C2H2/N20  could  also  be  approximated  for  conditions  beyond  the 
range  of  experimental  measurement.  This  correlation  will  be  used 
as  input  to  a  model  in  the  following  section  to  predict  theoretically 
another  dynamic  parameter  -  the  critical  initiation  energy. 

3.2.  C2H2/N20/Ar  critical  initiation  energy  measurement  and 
prediction 

Figure  6  shows  the  critical  energy  as  a  function  of  initial  condi¬ 
tions  obtained  from  both  experimental  measurement  and  theoret¬ 
ical  prediction  for  the  stoichiometric  C2H2/N20  mixture.  Similarly, 
the  upper  and  lower  limits  of  the  error  bar  represent  respectively 
the  last  energy  values  used  in  an  experiment  with  the  same  condi¬ 
tions  at  which  successful  and  failure  of  direct  initiation  are  ob¬ 
served.  This  representation  was  also  accepted  in  our  previous 
papers  on  the  problem  of  direct  initiation  of  detonation  [25,27]. 
In  general,  similar  to  common  fuel-oxygen  mixtures,  by  increasing 
the  initial  pressure  the  critical  energy  decreases.  For  the  theoretical 
prediction,  result  is  obtained  using  a  work  done  model  and  the 
dc  =  A  •  A i  critical  tube  correlation  obtained  in  Section  3.1.  The  work 
done  concept  [18,28,50,51]  is  essentially  based  on  a  piston  model. 
Assuming  the  planar  detonation  wave  as  an  interface  or  a  piston 
propagating  at  constant  velocity  and  transmitting  into  an  uncon¬ 
fined  space,  its  work  done  can  be  given  by: 

Ec=  [  pCjS(t)Uqdt  (2) 

Jo 

where  pc j,  Ucj  denote  the  CJ  detonation  pressure  and  particle  veloc¬ 
ity,  respectively,  t*  is  the  time  when  the  rarefaction  wave  reaches 
the  tube  axis,  t*  =  dc/2aQ  with  aQ  being  the  sound  speed  of  the  det¬ 
onation  products  [18,52].  S(t)  is  the  surface  area  of  the  planar  CJ 
detonation  wave  in  the  confined  tube,  which  is  written  as: 

S(t)  =  7r(y)2  (3) 


Equivalence  ratio  (9) 

Fig.  8.  Induction  length  as  a  function  of  equivalence  ratio  for  C2H2/N20  mixture  at 
p0  =  100  kPa. 


where  dc  is  the  tube  diameter.  Hence,  this  simplified  work  done 
model  gives: 


c  _Pquqftdc 
c  8aq 


(4) 


It  is  thus  suggested  by  Eq.  (4)  that  the  critical  tube  diameter  and 
critical  initiation  energy  can  be  directly  linked  with  each  other  by 
this  work  done  model.  With  the  knowledge  of  experimentally  mea¬ 
sured  and  theoretically  predicted  critical  tube  diameter  data  intro¬ 
duced  in  the  previous  section,  critical  initiation  energy  can  be 
approximated  and  compared  with  the  experimental  measurement. 
From  the  comparison  between  theoretical  prediction  and  experi¬ 
mentally  measured  data  which  is  shown  in  Fig.  6,  it  is  found  that 
the  theoretical  prediction  can  estimate  the  critical  energy  at  each 
initial  pressure  with  very  reasonable  accuracy. 

Figure  7  shows  the  variation  of  experimental  and  predicted  crit¬ 
ical  energies  versus  equivalence  ratio  for  C2H2/N20  mixtures  at  the 
initial  pressure  of  100  kPa.  Both  the  experimental  data  and  the  the¬ 
oretical  prediction  curves  indicate  that  the  critical  initiation  energy 
along  with  its  equivalence  ratio  cp  is  in  the  form  of  typical  “U” 
shape,  and  the  minimum  of  the  critical  initiation  energy  appears 
around  cp  =  2.5.  This  behavior  is  in  good  agreement  with  the 
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induction  zone  length  variation  curve  obtained  from  chemical 
kinetics  calculation,  which  is  shown  in  Fig.  8.  The  explosive  sensi¬ 
tivity  characterized  by  the  direct  initiation  energy  also  agrees  with 
the  critical  tube  diameter  results,  which  both  show  that  C2H2/N20 
mixture  is  most  sensitive  around  (p  =  2.5. 

To  look  at  the  dilution  effect,  Fig.  9  shows  the  critical  energy  for 
C2H2/N20  mixtures  diluted  with  different  amounts  of  inert  argon. 
As  seen  from  Fig.  9  that  with  an  increase  of  the  amount  of  argon 
dilution,  the  critical  initiation  energy  increases  consequently;  In 
other  words,  with  more  argon  dilution,  the  mixtures  tend  to  be  less 
sensitive  to  detonation.  As  shown  in  Fig.  10  the  variation  of  induc¬ 
tion  length  as  a  function  of  amount  of  Ar  dilution,  by  increasing  the 
amount  of  argon  dilution,  the  induction  zone  length  increases, 
which  can  also  explain  qualitatively  the  increase  of  critical  initia¬ 
tion  energy. 

3.3.  Critical  energy  comparison  between  C2H2/N20/Ar  and  H2/N20/Ar 
mixtures 

The  critical  energies  of  direct  detonation  initiation  for  stoichi¬ 
ometric  C2H2/N20  and  H2/N20  mixtures  as  a  function  of  initial 
pressure  are  compared  in  Fig.  11.  By  comparing  the  results  for 


Fig.  9.  Critical  energy  as  a  function  of  argon  dilution  from  experimental  measure¬ 
ment  and  theoretical  prediction  for  stoichiometric  C2H2/N20  mixture  at 
p0  =  1 00  kPa. 


Fig.  10.  Induction  length  scale  as  a  function  of  argon  dilution  for  stoichiometric 
C2H2/N20  mixture  at  p0  =  100  kPa. 


stoichiometric  C2H2/N20  and  H2/N20  mixtures  at  the  same  initial 
pressure,  it  is  clear  that  the  critical  energies  for  H2/N20  are  always 
higher  than  those  of  C2H2/N20  at  the  range  of  initial  pressures  con¬ 
sidered  in  this  study.  Figure  12  shows  the  results  of  induction 
length  scale  analysis;  As  can  be  seen  that  the  induction  length  of 
H2/N20  is  around  four  times  larger  than  that  of  C2H2/N20  mixture. 
Equivalently,  the  larger  induction  length  scale  results  in  the  bigger 
critical  energy  required  for  direct  initiation  of  a  detonation. 

The  experimentally  measured  critical  initiation  energies  for 
C2H2/N20  and  H2/N20  mixtures  at  atmospheric  pressure  with  the 
variation  of  equivalence  ratio  are  shown  in  Fig.  13.  Both  curves 
of  critical  energy  variation  with  equivalence  ratio  for  C2H2/N20 
and  H2/N20  mixtures  are  again  described  by  the  typical  “U”  shape 
behavior.  However,  the  minimum  critical  energy  of  H2/N20  is 
around  cp  =  0.8,  and  for  the  C2H2/N20  mixtures,  the  minimum  crit¬ 
ical  energy  occurs  around  (p  =  2.5.  It  also  appears  that  the  C2H2/N20 
is  more  sensitive  than  H2/N20  at  the  rich  side.  This  scenario  is  very 
similar  from  the  chemical  kinetic  point-of-view  demonstrated  by 
the  induction  length  variation  with  equivalence  ratio  for  both  mix¬ 
tures,  which  are  shown  in  Figs.  8  and  14.  It  can  be  seen  that  the 
minimum  critical  energy  always  occurs  at  the  composition  where 
induction  length  has  the  minimum  value. 


Fig.  11.  Critical  energy  as  a  function  of  initial  pressure  for  stoichiometric  C2H2/N20 
and  H2/N20  mixtures. 


Fig.  12.  Induction  length  as  a  function  of  initial  pressure  for  stoichiometric  C2H2/ 
N20  and  H2/N20  mixtures. 
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Fig.  13.  Critical  energy  as  a  function  of  equivalence  ratio  for  (a)  H2/N20  and  (b)  C2H2/N20  mixtures  at  p0  =  100  kPa. 
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Fig.  14.  Induction  length  as  a  function  of  equivalence  ratio  for  H2/N20  mixture  at 
p0=  100  kPa. 
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Fig.  15.  Critical  energy  as  a  function  of  %  argon  dilution  for  stoichiometric  H2/N20 
and  C2H2/N20  mixtures  at  the  initial  pressure  of  p0  =  100  kPa. 


The  critical  initiation  energy  results  for  stoichiometric  C2H 2\ 
N20  and  H2/N20  with  different  amounts  of  argon  dilution  mixtures 
are  shown  in  Fig.  15.  One  can  note  that  with  an  increase  of  argon 
dilution,  the  critical  energy  for  both  mixtures  increases.  However, 
it  increases  only  very  slightly  for  C2H2/N20/Ar  mixture,  yet  for 
H2/N20/Ar  mixture  the  increase  of  critical  energy  appears  much 
steeper.  In  addition,  the  results  of  induction  length  analysis  are 
shown  in  Fig.  16,  which  indicate  that  for  the  H2/N20  mixture,  the 
induction  zone  length  scale  is  always  bigger  than  C2H2/N20  when 
diluted  with  the  same  amount  of  argon. 

One  interesting  observation  is  that  although  the  ZND  induction 
lengths  for  both  mixtures  do  not  vary  significantly  with  increasing 
argon  dilution,  one  can  notice  in  Fig.  15  that  there  appears  a  fast 
increase  in  critical  energy  for  H2/N20  mixtures  with  dilution  above 
30%  and  for  C2H2/N20  mixtures  diluted  with  more  than  50%.  This 
fast  increase  in  the  critical  energy  could  perhaps  be  related  to 
the  dynamic  effect  instead  of  the  chemical  effect  [25,26,53].  It  is 
known  that  by  increasing  the  argon  dilution,  the  detonation  wave 
becomes  more  stable  in  the  sense  that  the  reaction  zone  structure 
is  more  regular  with  weaker  transverse  waves  [54],  and  the  dy¬ 
namic  change  in  the  regularity  of  the  detonation  structure  is 
known  to  have  an  effect  on  the  propagation,  limits  and  initiation 


of  a  detonation  wave  [55].  Previous  studies  proposed  that  the  sta¬ 
bility  of  the  detonation  wave  can  be  characterized  by  a  stability 
parameter  calculated  by  the  effective  activation  energy  of  the 
induction  zone  multiplied  by  the  ratio  of  induction  to  exothermic 
heat  release  length  [53,56].  As  shown  in  our  previous  work  [25], 
stable  or  regular  mixtures  characterized  by  small  instability 
parameter  %  tend  to  be  less  sensitive  to  detonation  initiation.  Look¬ 
ing  at  Fig.  17  which  shows  the  stability  parameter  %  for  both  mix¬ 
tures  as  a  function  of  the  degree  of  argon  dilution,  by  increasing  the 
argon  dilution  the  stability  parameters  x  f°r  both  mixtures  de¬ 
crease  and  the  detonation  shall  become  more  stable  and  therefore, 
the  increase  in  critical  energy  as  shown  in  Fig.  15  can  be  explained 
from  this  detonation  instability  consideration.  Furthermore,  the 
rate  of  decrease  in  x  for  H2/N20  mixtures  by  increasing  the  argon 
dilution  appears  to  be  faster  than  that  of  C2H2/N20,  and  hence,  it 
also  explains  why  the  onset  of  the  fast  increase  in  critical  energy 
for  H2/N20  mixtures  is  earlier  at  about  30%. 

From  the  critical  energy  comparison  between  C2H2/N20/Ar  and 
H2/N20/Ar  mixtures,  all  results  suggest  that  C2H2/N20/Ar  is  more 
detonation  sensitive  than  H2/N20/Ar  at  the  same  initial  conditions, 
which  can  also  be  qualitatively  supported  by  the  ZND  analysis  that 
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■  C2H2-N20  experimental  data 

•  C2H2-02  experimental  data 

- C  H  -N  O  experimental  data  curve  fit 


Fig.  16.  Induction  length  as  a  function  of  %  argon  dilution  for  stoichiometric  H2/N20 
and  C2H2/N20  mixtures  at  p0  =  100  kPa. 


Ar,  (%) 

Fig.  17.  Stability  parameter  x  as  a  function  of  the  degree  of  argon  dilution  for 
stoichiometric  H2/N20  and  C2H2/N20  mixtures. 


the  induction  length  scale  for  C2H2/N20/Ar  is  always  smaller  than 
that  of  H2/N20/Ar  mixture. 


3.4.  Critical  energy  comparison  between  C2H2/N20/Ar  and  C2H2/02/Ar 
mixtures 

To  study  the  effect  of  the  oxidizer  on  the  critical  initiation  en¬ 
ergy,  a  series  of  experiments  using  mixtures  of  C2H2/02/Ar  at  the 
similar  initial  conditions  as  those  with  C2H2/N20/Ar  mixtures  were 
also  performed  in  order  to  make  comparison  of  their  critical  ener¬ 
gies  and  further  investigate  chemical  kinetic  effects  in  using  differ¬ 
ent  oxidants. 

Figure  18  shows  the  critical  energy  of  direct  detonation  initia¬ 
tion  for  stoichiometric  C2H2/N20  and  C2H2/02  mixtures  as  a  func¬ 
tion  of  the  initial  pressure.  Due  to  the  different  order  of 
magnitude  of  the  detonation  sensitivity  of  these  mixtures,  and  to 
guarantee  the  accuracy  of  the  measurement,  the  critical  energies 
of  both  mixtures  were  measured  at  the  different  range  of  initial 
pressure.  Nevertheless,  from  the  tendency  of  the  curve  fit  with 
experimental  data,  it  is  readily  seen  that  the  critical  energy  C2H 2\ 
N20  mixture  is  around  2  orders  of  magnitude  bigger  than  that  of 
C2H2/02  if  the  same  initial  pressure  is  considered.  Once  again,  from 


Initial  Pressure  (kPa) 

Fig.  18.  Critical  energy  as  a  function  of  initial  pressure  for  stoichiometric  C2H2/N20 
and  C2H2/02  mixtures. 


Fig.  19.  Induction  length  as  a  function  of  initial  pressure  for  stoichiometric  C2H2/ 
N20  and  C2H2/02  mixtures. 


the  aspect  of  induction  length  analysis,  which  is  shown  in  Fig.  19,  it 
is  observed  that  the  ZND  induction  length  scale  for  C2H2/N20  is 
around  three  times  bigger  than  that  of  C2H2/02  when  at  the  low 
pressure,  and  about  the  ratio  is  about  10  as  the  pressure  increases 
above  100  kPa.  It  is  also  interesting  to  notice  that  the  ratio  between 
the  two  critical  energies  of  C2H2/N20  and  C2H2/02  mixtures  found 
in  Fig.  18  is  of  the  order  of  400.  Using  the  ZND  induction  length  as 
the  characteristic  length  scale,  from  dimensional  consideration  the 
initiation  energy  should  scale  according  to  Ec  ~  A3  and  hence,  the 
ratio  of  the  ZND  induction  length  between  these  two  mixtures 
should  be  approximately  of  the  order  of  7.  This  agrees  well  with 
the  result  shown  in  Fig.  19,  which  shows  on  average  a  ratio  of  6 
over  the  initial  pressure  range  considered  in  this  study.  A  discrep¬ 
ancy  is  expected  since  it  is  known  that  the  scaling  Ec  ~  A3  is  not 
appropriate  for  unstable  detonations  with  irregular  reaction  zone 
structures  and  the  energy  formula  in  general  depends  also  on  the 
CJ  detonation  velocity,  density  p0,  specific  heat  ratio  y,  cell  size  2, 
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■  C2H2-N20  experimental  data(p0=100kPa) 


Fig.  20.  Critical  energy  as  a  function  of  equivalence  ratio  for  C2H2/N20  and  C2H2/02 
mixtures. 


Fig.  21.  Induction  length  as  a  function  of  equivalence  ratio  for  C2H2/N20  and  C2H2/ 
02  mixtures  at  p0  =  1 00  kPa. 


etc.,  Ec  ~  PoDcj^j  [57].  These  parameters  vary  also  between  C2H 2\ 
N20/Ar  and  C2H2/02/Ar  mixtures. 

The  critical  initiation  energies  for  C2H2/N20  and  C2H2/02  mix¬ 
tures  at  different  equivalence  ratios  were  also  experimentally  mea¬ 
sured  to  observe  whether  the  minima  move  by  using  different 
oxidizers  and  results  are  shown  in  Fig.  20.  It  should  be  noted  that 
the  initial  pressure  does  not  change  the  minimum  critical  energy 
behavior  with  the  variation  of  equivalence  ratio,  and  again  due  to 
the  high  detonation  sensitivity  of  C2H2/02  mixtures  and  experi¬ 
mental  limitations,  the  initial  pressure  of  C2H2/02  of  10  kPa  is  used 
instead  in  this  study.  It  can  be  seen  from  Fig.  21  that  the  minimum 
critical  energies  locate  at  about  cp  =  2.5  for  both  mixtures.  If  one 
looks  at  the  induction  zone  length  variation  with  equivalence  ratio, 
the  induction  zone  length  also  arrives  at  the  minimum  value 
around  (p  =  2.5  for  both  mixtures  under  the  atmospheric  pressure. 
The  result  shows  that  the  minimum  critical  energy  behavior  is  in 
agreement  with  that  of  induction  zone  length. 

The  last  part  of  this  study  aims  to  look  at  the  critical  initiation 
energy  for  stoichiometric  C2H2/N20  and  C2H2/02  mixtures  diluted 
with  different  amounts  of  argon  and  the  results  are  shown  in 
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Fig.  22.  Critical  energy  as  a  function  of  %  argon  dilution  for  stoichiometric  C2H2/N20 
and  C2H2/02  mixtures. 


Fig.  23.  Induction  length  as  a  function  of  %  argon  dilution  for  stoichiometric  C2H2/ 
N20  and  C2H2/02  mixtures. 


Fig.  22.  For  the  undiluted  case  and  22%  argon  dilution  in  C2H2/02 
mixtures,  the  critical  energies  were  obtained  by  the  theoretical 
predictions  which  are  based  on  the  experimental  measured  cell 
sizes  data  combined  with  the  model  presented  earlier  in  this  paper 
and  the  classical  relationship  dc  =  132.  The  cell  size  for  the  undi¬ 
luted  and  22%  argon  diluted  C2H2/02  mixtures  are  extracted  from 
the  CALTECH  detonation  database  [58]  and  the  dissertation  by  Rad- 
ulescu  [53],  respectively.  The  experimental  data  curve  fits  show 
that  with  increasing  amount  of  argon  dilution,  the  critical  energies 
increase  for  both  mixtures.  However,  the  critical  energy  for  C2H 2\ 
N20/Ar  mixture  is  always  2  orders  of  magnitude  larger  than  that 
of  C2H2/02/Ar  mixture  with  the  same  amount  of  argon  dilution. 
Figure  23  shows  the  induction  length  variation  for  these  mixtures; 
it  is  apparent  that  the  induction  length  for  C2H2/N20/Ar  is  around 
one  order  of  magnitude  bigger  than  that  of  C2H2/02/Ar  mixture. 

From  the  critical  energy  comparison  between  C2H2/N20/Ar  and 
C2H2/02/Ar  mixtures,  it  can  be  found  that  the  mixture  tends  to  be 
less  detonation  sensitive  by  using  nitrous  oxide  as  oxidizer  than 
using  oxygen,  thus  more  energy  is  needed  to  achieve  a  successful 
direct  initiation.  The  experiment  results  also  indicate  that  the  oxi¬ 
dizer  does  not  change  the  overall  critical  energy  behavior  and  the 
location  of  its  minimum. 
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4.  Concluding  remarks 

In  this  study,  detonation  critical  tube  diameters  were  measured 
for  C2H2/N20/Ar  mixtures  at  various  initial  conditions  (i.e., 
(p  =  0.625-2.5  and  argon  dilution  up  to  50%,  respectively).  Using 
chemical  kinetic  calculation  with  Konnov  mechanism  (Release 
0.3)  and  present  experimental  data,  a  simple  correlation  for  the 
detonation  critical  tube  diameter  of  C2H2/N20/Ar  mixtures  with 
the  ZND  induction  zone  length  scale,  i.e.,  dc  =  594.8<p0623 
(1  -XAr)02176(p/p°)“00246d/  is  developed  valid  for  the  stoichiome¬ 
tric  mixture  at  different  initial  pressures  ranging  from  50  to 
130  kPa,  equivalence  ratios  from  0.625  to  2.5  and  stoichiometric 
mixtures  with  the  maximum  percentage  of  Ar  dilution  up  to  50% 
at  atmospheric  pressure.  The  critical  energy  for  direct  blast  initia¬ 
tion  can  also  be  predicted  by  combining  this  correlation  function 
with  a  work  done  model,  and  the  results  show  reasonable  accuracy 
by  the  comparison  with  those  measured  experimentally. 

The  critical  energy  results  for  direct  initiation  of  detonation  in 
C2H2/N20/Ar  mixtures  are  compared  with  those  in  H2/N20/Ar  mix¬ 
tures.  The  results  indicate  C2H2/N20/Ar  is  more  sensitive  than  H 2\ 
N20/Ar  at  the  same  initial  conditions.  The  effect  of  the  oxidizer 
on  the  detonation  sensitivity  is  also  investigated  by  comparing 
the  critical  energy  between  C2H2/N20/Ar  and  C2H2/02/Ar  mixtures. 
The  results  show  that  the  mixture  tends  to  be  less  detonation  sen¬ 
sitive  by  using  nitrous  oxide  as  oxidizer  than  oxygen  if  the  same 
initial  condition  is  considered,  thus  higher  initiation  energy  is  re¬ 
quired  to  achieve  direct  initiation  of  a  spherical  detonation  in 
C2H2/N20/Ar  than  in  C2H2/02/Ar  mixtures.  To  explain  qualitatively 
these  variations  of  critical  initiation  energy  and  critical  tube  diam¬ 
eter  observed  in  the  experiments,  ZND  induction  length  analysis  is 
also  performed.  Such  analysis,  from  chemical  kinetic  consideration, 
has  proven  to  provide  some  insights  on  the  qualitative  trends  of 
detonation  dynamic  parameters  with  different  initial  conditions. 
In  all  cases,  the  critical  initiation  energy  behavior  agrees  accord¬ 
ingly  with  the  variation  of  the  ZND  induction  length. 
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